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Abstract 



Given a real valued function on R™ we study the problem of recovering the function from 
its spherical means over spheres centered on a hyperplane. An old paper of Bukhgeim and 
Kardakov derived an inversion formula for the odd n case with great simplicity and economy. 
We apply their method to derive an inversion formula for the even n case. A feature of our 
inversion formula, for the even n case, is that it does not require the Fourier transform of the 
mean values or the use of the Hilbert transform, unlike the previously known inversion formulas 
for the even n case. Along the way, we extend the isometry identity of Bukhgeim and Kardakov 
for odd n, for solutions of the wave equation, to the even n case. 



1 Introduction 



Below n > 1 will be a positive integer, a point x 6 M. n will occasionally be written 

with x' 6 M n , x n € M, and \S n I = — -. — j-^r will be the surface area of the unit sphere in M". We 

r(n/2) 

take /(£) = L n e~ lx '^f(x) dx as the definition of the Fourier transform of /. Define C^°(M n ) to be 



the set of smooth functions f(x) on R n which vanis fill to infinite order on x n = and are even in 
x n . Let Cg^(R n ) denote the subspace of compactly supported functions in C^°(R n ). Finally, define 
L 2 , w (W l ) to be the Hilbert space of measurable functions f{x) on R n which are even in x n and for 

which J Rn dx is finite - we use the associated inner product. 

For any continuous function f(x) on R n define the spherical mean value operator 

(Mf)(x,t) = — ^ / f(x + t0)d9, (x,i)GR n xR; (1.1) 
\S I J\e\=i 

(M f)(x,t) is the spherical average of / on the sphere centered at x, radius t. The goal is the 
recovery of /, given (M f)(x,t) for all x on the hyperplane x n = and for all t G R. Clearly on 
x n = 0, (Mf)(x,t) does not change if f{x) is replaced by f(x', —x n ), its reflection across x n = 0. 
Hence one can recover at most the even, in x n , part of /. If / G C^(R n ) then (M f)(x,t) is an 
even function of x n and t and (M f)(x' ,0,t) vanishes to infinite order on t = 0. So we define the 
map M : C^(M. n ) q?°(R n ) by 

CM/)(x', t) := (M/)(x', 0, t), (x', t) G 

The goal of this article is to construct an inverse of Ai and characterize the range of M. We will 
do this indirectly. 

There is a well known relationship between spherical mean values and solutions of the wave 
equation. For any / G C^(W l ), let u(x,t) be the solution of the initial value problem 

u tt - Au = 0, on 1" x R, (1.2) 
u (x,0) = f(x), u t (x,0) = 0, x G R n . (1.3) 

Note that u(x,t) is even in i and also in x n because f{x) is even in x n . From page 682 of [CH62] 
there is a formula^ for u(x, t) in terms of Mf and in terms of /. For any smooth even function h(t) 
on R define the differential operator 

(DW) := If, t e R. 

Note that Dh(t) makes sense for all t because h is even and {Dh){t) is an even function of t. 
Further, D may be considered as differentiation with respect to t 2 because it may be verified that 
D(h(t 2 )) = h'{t 2 ). We have the following relations for all (x,t) G R n+1 , t / 0: for odd n 

U{x ' t] = f^2)t D(n ~ 1)/2 ( tn ~ 2 ( M f)^t)) (I-*) 
= 7r (l-n)/2 tD (n-i)/2 f 6 (t 2 -\x-y\ 2 )f(y)dy, (1.5) 

1 that is (d a f)(x) = on x n = for all multi-indices a 
there seems to be a misprint in the formula in [CH62] for the even n case 
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and for even n 

«om) = fm tDnl2 Slw=^ Mf) ^ r)dr (L6) 

= ^HD^( H f~\*~ y ^ f(y)dy (1.7) 

where -ff(s) is the Heaviside function. Using (jl.2j) and that / vanishes to infinite order on x n = 0, 
it may be verified that u(x' , 0, t) vanishes to infinite order on t = 0. Define the map U : C^(M. n ) — > 
C™(R n ) by 

(lAf){x',t) :=u(x',0,t). 

We construct the inverse of IA and characterize the range of IA which will give us the inverse of A4 
and implicitly characterize the range of Ai. 

Define the tempered distribution K(t) := t~ 1 / 2 H(t), t G R and let K^ m \t) denote the m-th 
derivative of K(t). 



Theorem 1. Suppose n > 1 is even. 

(a) (Isometry) For any f G Co^(M n ) we have the isometry 



W'fcfipfrt, (i, 



\ x n\ JR n \t, 

solA extends to an isometry on L 2 w (MJ l ). 

(b) (Inversion Formula) For any f G Cq^(W 1 ), x G M. n , x n > we have 

(c) (Range Characterization) The range of IA is a closed subset of L 2 w (W l ). Further, if F G 
L 2 w (W n ) is in the range of the restricted IA (domain is C^iW 1 )) then 

(i) F(x',t) G C™(R n ) and H(t)F(x', y/r)/y/r G H s (W n ) for all real s, 

(ii) \d%d{F(x',t)\ < C j>a (l + \t\)-^ n ^/ 2 for all j > and all multi-indices a. 

Finally, if F satisfies (i), (ii) then F is in the range ofU iff'ix G R n , x n ^ 0, we have 

H(* + *l-\yf) t (L^^^t) ^ 



JR™- 1 J 



y/t 2 + xl - \y'\ 2 \2tdt 



In the statement of Theorem [T] we were careful to write only absolutely convergent integrals 
which forced us into a long-winded statement of (c). A briefer version of (b) and (c) would be that 
for all / G Ll w (R n ) 

f( x ) = ^k[ K^ 2 \t 2 -x 2 n -\y'\ 2 )(Uf)(x' + y',t)dy'dt, \fx G R n , x n > 0, (1.11) 
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and an F G 1^(1") is in the range of U iff 

f K^ 2 \t 2 + x 2 n -\y'\ 2 )F{x' + y' ,t)dy' dt = 0, Vx G R n , x n / 0. (1.12) 

JR™ 

However, now (II. lip and (jl . 12H have to be interpreted in a sense buried in the proof of Theorem [TJ 
The result for odd n given in [BK78] has a form similar to the formal expressions above. They can 
be rewritten in terms of absolutely convergent integrals similar to the expressions in Theorem [TJ 

From the hypothesis of Theorem [TJ it would seem that the inversion formula is not applicable to 
an / G C^°(R n ) which is even in x n but which does not have the right decay near x n = 0. However, 
for such /, if u(x',0,t) is known then using u Xn (x', 0, i) = (because u is even in x n ) and f)1.2|) 
one may determine /(x',0) = u(x',0,0) and all the x n derivatives of f{x) = u(x,0) on x n = 0. 
Knowing these derivatives of f(x) on x n = we can create a known function g{x) G C^°(lR n ) so 
that f = f — g vanishes to a chosen finite order (or infinite order) on x n = 0. Since g is known 
we can solve (jl.2p . (jl.3p with g replacing / and hence from u(x',0,t) we can obtain the trace on 
x n = of the solution of (|1.2p . (|1.3p with / replaced by /. Now we can apply the inversion formula 
to the / data and recover / and hence /. 

The early work on the inversion of IA (and hence of A4) is described in [Jo55]. For odd n, a 
nice inversion formula for IA (and hence M) was given by Bukhgeim and Kardakov [BK78]. Later, 
inversion formulas for M for odd and even n were derived in [No80], [Fa85], [An88], [NRT95], 
[K103], and [Be09]; [SQ05] has a numerical inversion scheme. The problem under consideration 
has applications to inverse problems in elasticity and other areas - see [BK78] and [Fa85]. The 
article [FPR04] studies the similar problem of recovering a compactly supported function from its 
spherical averages over spheres centered on the boundary of a region containing the support of the 
function. 

Beltukov's formula in [Be09], for the odd n case, is essentially the inversion formula in [BK78]. 
Beltukov rediscovered some of the ideas in [BK78] (he seems not to have been aware of [BK78]) 
but not all the ideas in [BK78]. Applying the ideas in [BK78] to the even n case, we give a concise 
derivation of an inversion formula for the even n case which does not require the Fourier transform 
of the data or the use of the Hilbert transform. In [BK78] an isometry identity was established for 
solutions of the wave equation for the odd n case; we extend that identity to the even n case. 

The proof of Theorem [TJ uses the tempered distribution N(x) := tt~ ti I 2 K^ n ^ 2 \x n — \x'\ 2 ) and 
its Fourier transform. 

Proposition 2. For even n > 1 we have N(£) = e l ^'l /(^™) for all £ G M. n . 

To prove (c) in Theorem [TJ and to be sure that the integral in (b) is absolutely convergent we 
need some decay estimates for solutions of the wave equation - see Theorem 1.1 in [So08]. 
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Proposition 3. Suppose n > 1 is a positive integer, f,g G C^°(M. n ) and v(x,t) the solution of 

v tt -Av = 0, on M. n x R, (1.13) 

v(x,0) = f(x), v t (x,0) = g(x), (1.14) 

then 

\v(x,t)\ 4 {l + \t\)- (n ^/ 2 , n odd 

\v(x,t)\ 4 (l + |t|)~ (n ~ 1)/2 (l + ||t| - M|)~ (n ~ 1)/2 , n even (1.15) 
with the constant determined by f,g,n. 

For n odd, v(x,t) is supported in \t — \x\ \ < R if /, g are supported in \x\ < R. This is not true 
for n even, instead we have an extra decay contribution over the region \t — \x\\ > R. We have 
included the statement and the proof for the odd n case because of its possible application to the 
inversion formula for odd n in [BK78]. 

The rest of the article is as follows: Section [2] contains the proof of Theorem[Tl Section[3]contains 
the proof of Proposition [21 and Section U] contains the proof of Proposition [3l We have included 
the proof of Proposition [3] for completeness though it is stated as a theorem in [So08] whose proof 
was left as a good exercise. The proof of Theorem [1] would be quite short (see [BK78] for the odd 
n case) if we accept formal arguments and expressions; our proof is somewhat long because of the 
need to show the absolute convergence of the integrals in the expressions. 

This work was done while the first named author was supported in part by a grant from 
the UGC via DSA-SAP and the second named author was on sabbatical at the Department of 
Mathematics of the Indian Institute of Science in Bangalore, India. The second named author 
thanks the department for its warm hospitality and generous financial support and the University 
of Delaware for granting the sabbatical. 



2 Proof of Theorem [T] 

For t > and x' £ R n_1 from (|1 .7j) we have (integrals below to be interpreted as distributional 
action) 

(Uf)(x',0,t) = U (x',0,t) =7T _ n/2 I K (n/2 ){t 2_ W _ yf _ y 2 n)f{y)dy 

t t J^n 

= 2VT-/ 2 / / K^\t 2 - \x> -yf- yl) f(y) dy n dy> 

JR"- 1 JO 

-nil f f°° TAnli\/j2 I / /|2 \ f(l^VOn) , , / 
= 7T 1 I I K K 1 >(t -\x - q\ - q n ) -£z — dq n dq . 

Jr"- 1 Jo \/Qn 
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Hence taking x' = p' , t = ^/p^ for any p' E W 1 , p n > 0, we have 

(^f)\Pi\/Pn) - n /2 f f°° T An/2)f i / /|2 \ /(# ' V^n) A i 
=- = vr ' / / K y 1 >(p n - \p - q | - q n ) — dq n dq . (2.1) 

JR™- 1 JO vQn 



Let L+(R n ) denote the subspace of functions g(x) in L 2 (R n ) which are zero for x n < 0. Define 
the specialized zero conversion operator Z : L 2 w (M. n ) — > L^_(R n ) with 

[0 p n < 0. 

Then i? is an isomorphism from Lg U) (]R n ) to L^(R n ) and its inverse is given by (Z~ 1 g)(x > ,x n ) = 
x n g(x' ,x n ). Further Z maps C y ^(W l ) to the subspace of functions in Co°(IR n ) which are zero on 

X n < 0. 

For any / 6 C^iW 1 ) let g = Zf; then from flZI]) we have 
(iJZ^)(p)=^/ 2 / ^^(Pn-lp'-g'l 2 -^)^)^^^*^^), PGM"; (2.2) 

note that this is true even if p n < 0. This suggests we define a map N : Co°(M n ) — > C°°(M n ) with 

(jVg)( P ) = (N*g)(p), P eR n - 

then -ZZY/ = MZf. Using Proposition [2] and the definition of M we have 

WO = HMO = ^ mn) 9{0- (2-3) 

So \ {Ng){C)\ = |<?(£)l f° r an £ ^ ^™ anc ^ N extends to an isometric isomorphism from L 2 (R n ) to 
L 2 (W n ). In particular, for all / G Cq^OR 71 ) we have ||^/||l 2 (r™) = ll^^/IU 2 (R") an d hence 

/ / aPn dp = / * dg n d<? . 

JR™- 1 JO Pn JR"- 1 JO Qn 

Introducing the change of variables p' = x',p n = x n , and y' = q',t = yjq^ on the regions p n > 
and q n > we obtain (|1.8p which proves (a) of Theorem [TJ 

The equation (12.3P allows us to express any g € L 2 (M n ) in terms of A/"<7 giving us 

m = e -i\i'\*Mn)jfijft = -S n )jfij($. (2.4) 

Now, for any g £ Co°(IR n ), g(£) is a rapidly decaying function, and N(£', £ n ) is bounded, so Ng(£) is 
a rapidly decaying function and hence Mg £ H m (M. n ) for all real m. Further N(£', —£, n ) is bounded 
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so N(q' , —q n ) G H m (M. n ) for any real m > n/2; pick one m. Then for any p G W 1 , using the 
Fourier inversion formula and (|2.4p . we have (below T p is the translation operator) 

(27r) n J Rn {2n) n J R n 

= 7T^ I Vv^(0^= {N(-J,q n ),(Mg)(p + q)) 

\ Z7r ) Jut™ 

= {N( J /-J,q n -p n ),(Ar g )(q)) (2.5) 
where the (, ) is to be understood as the action of an element in H~ m (M n ) on an element of H m (M n ). 

For any / G C^(M n ), take g = Z/, then A/"g = ZUf; hence for any p G M n , from (|2.5p we have 

(Zf)(p) = (N(jt - q', q n - Pn ), (ZUf)(q)). (2.6) 

Working formally, for the moment, (12. 6h implies 

(Zf)(p) = (K^ 2 \p'-q',q n -p n ),(ZUf)(q)) 

/ a \ n/2 

= (-1)"' ! {A>' - ,',,„ -p„), ( ^-J (£«/)(«)) 

JM n \fq n - Pn~ W - P\ \ a 1-n./ 

For any i£l™ with x n > 0, taking p! = x',p n = x\, and taking y' = q',t = ^fq~^ (note g n > on 
the support of ZUf), we have 

/(*)_o, _W2 /" r ^(t 2 -^-|x'-yf)^l 5^"/ 2 (^/)(y',t) 



thus proving (b) of Theorem [TJ 

We now address (c) of Theorem [H If / G Cf£(R n ) then (Uf){x',t) G C~(R n ) and from ^ 
G -£P(R n ) for all real s. Further, if u is the solution of (|1.2|) . (|1.3p then <9"<9 f m u is a solution 
of (|1.13p . (|1.14p for some f,g G C£°(R n ) and hence from Proposition [3] 

i^ar(^/)(x',t)i < c m>a (i + |tir (n ~ 1)/2 

for all m > and all multi-indices a, proving (i), (ii) of (c) in Theorem [TJ 

The extended U is an isometry and hence its range is closed. A function F(x',t) G Ll w (M. n ) 
is in the range of U iff F = Uf for some / G L 2 w {W l ), that is (using the isomorphism of Z) iff 
ZF = ZUf = ATZf, that is (using the isomorphism of AT) iff Zf = N~ X ZF for some / G L^ w (R n ). 
So a necessary condition for F to be in the range of U is that (Af- 1 ZF)(p) = for p n < 0. 
Conversely, if (J\f~ 1 ZF)(p) = on the region p n < 0, then g := N~ 1 ZF is in L\{R n ). Let 
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/ := Z~ l g- then / £ h\ w (R n ) and 3 = Zf and hence Z/ = N~ X ZF implying ZF = MZf = ZUf, 
so F = Uf proving our claim. 



If F satisfies (i), (ii) of (c) in Theorem[H then ZF £ H s (R n ) for all s. From fl23D, (N~ 1 ZF)(£) = 
N(£' , —^ n )ZF(^) and since the right side is an integrable function so is the left side, and hence by 
the Fourier inversion formula, as done for l\2.5\) . we have 

(J\f~ l ZF)(p) = (N(p' - q',q n - Pn ), (ZF)(q)), p £ R n . 

Hence, working formally for the moment, F is in the range only if for all p' £ M n_1 , with p n > 
we have (as done earlier) 

= {Af~ 1 ZF)(p , —pn) = (N(p' -q',q n +p n ),(ZF)(q)) 

1 \nl2 [ H(q n +p n -\q' -p'\ 2 ) ( d \ n/2 , 
= ( _7r ) / / 1 , ,19 ^~ {ZF)(q,q n )dq 

J^ n V1n+Pn-\q-P\ 

_ r r m^Kzm tff mb&> ^ M . 

Jr«-i Jo Vg„ + p n - W - q'\ 2 \oq n J v 7 ^ 

For any x £ K n with x n 7^ 0, take p' = x' , p n = x\. Also, inside the integral use the change of 
variables y' = q' , t 2 = q n . Then F is in the range of U iff for all x £ W 1 with x n 7^ we have 

Jk—Jo s/t2 + xl-\x'-y'\'' \2tdt) t 
proving (c) of Theorem [TJ 



So to complete a rigorous proof of the Theorem, for both parts (b) and (c), we have to show 
that if F satisfies (i), (ii) of part (c) of Theorem [TJ then for all p £ M. n , we have 

(N(p>-q',q n -p n ),(ZF)(q)) = (-vr)™/ 2 / H ^ ~ Pn ~ W?) M V^^^y^j^ 

JR n \Jqn-Pn-\q\ \°9n/ 

where the inner product (-, •) is to be interpreted as the action of an element in H~ m (W l ) on an 
element of H m (W n ). We now give a proof of this claim. 



For any <j> £ C yo (W n ), the value of (N(p'—q', q n —p n ), <t>{<l))i as determined by the duality between 
H~ m and H m , is the same as the value of the action of the tempered distribution N(p' — q', q n — p n ) 
on <j){q). Let x G C°°(R) with x(s) = 1 for s < and xO) = for s > 1 and if) £ C 00 ^"" 1 ) with 
ip(x') = 1 if \x'\ < 1 and ip(x') = if \x'\ > 2. Then 

MQ) ■= x(log(^) - k) ifj(q' /k)(ZF)(q) 
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X(s) 




Figure 1: Graphs of cft(s) and tp(x') 



is in C^°(M n ) and <f>k{q) converges to (ZF)(q) in H m (M. n ) for all m. Hence 
(N(p' - q,q n -Pn), {ZF){q)) = lim (N(p' - q' ', q n - p n ) , <t> k (q)) 



k— ¥00 



lim {-K) n l 2 {K{q n -p n - \p' - qf), {d/dq n ) n / 2 4> k {q)) 

k— >oo 

lim (-Tr)"^ ( n/2 )(K(q n ~Pn ~ \p'-q'\ 2 ),di( X (log(q 2 n )-k))^q'/k)d^(ZF)(q)) 
j=0 V J 7 

lim {-^{ K {q n -p n - \p' - q'\% X (log(q 2 n ) - k) ^(q' /k) d?/ 2 (ZF)(q)) (2.7) 

k— >oo 



m/2 f H(q n -p n -\q\ 2 ) ( d \ 1 , , 

-k) ' I , ==- -7T- (ZF)(q +p ,q n )dq. 

\fqn-Pn-Wr \ a< lnJ 



(2. 



The equations above required the observation that 



lim (K(q n -p n - \p' - q'\ 2 ), ^(x(log(^) - k)) ^{q> /k) d^(ZF)(q)) =0, j > (2.9) 

k— >oo 



and that (|2.7p equals (|2.8p . Both these observations follow from the use of the dominated conver- 
gence theorem, that for j > 0, d 3 n {xi^og{q 2 l ) — k)) is supported in e fc//2 < \q n \ < e^ k+1 ^ 2 , and the 
integrability in q of K{q n — p n — \p' — q'\ 2 ) (1 + \q n \)~^ d^ 2 '^ (ZF)(q) for j = • • • re/2. 



We now prove the integrability assertion in the previous sentence. For any function h{t) and 

III 

2tdt [ 



any real number k we have ^f'm{t k h{t)) = c\t k 2 h{t) + C2t k 1 dth(t), so for any m > we have 



1 1 \ m m m 

Ttdt ) {tkk{t)) = S cit k ~ l ~ 2{m ~ l) d l Mt) = E cit k ' 2m+l d\h{t). (2.10) 



Hence using the substitution q n = t 2 we have 
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and from the decay properties of F we conclude that 



m 



£i ww 



=( (1 + Itl)- 1 -"'!! + |i|)" l "" 1,/2 | ( .^ * (l + l9„l)- (2m+ " +1)/ ''. (2.11) 



Hence for any j > 0, 

7 „'|2<i 



/ / , | = jf (l + kn|) 3 dg d?. 



Pn ^kT<gn-pn V <fo ~ Pn ~ W \ 2 

z' + v' . a 

(l + Wn\r j (q n -Pn) in - 2)/2 dz'dq n 



x , i , d n n/2) ~ j (ZF) ) ( Vq^K z'+p', g, 




Pn J\z'\ 2 <l ^/T 



\z'\ 2 



Pn J\z'\ 2 <l V 1 - \ Z \ 



>V 



which is finite; we have used (|2.1ip and noted that for all j = • • • n/2 

{-j + (n - 2)/2) - (n + 1 + 2(n/2 - j))/4 = -j/2 - 5/4 < -5/4. 



QED 



3 Proof of Proposition [2] 



Let 5(R n ) denote the usual space of smooth, rapidly decreasing functions on W 1 . Choose x( s ) 6 
C°°(R) with x( s ) = e ~ s f° r s > and x( s ) = f° r s < — lj then % e , defined through 

satisfies lim e ^ + x e * <A = for all G 5(M n ) in the topology of 5(M n ). 



If P(x) = K(x n - \x'\ 2 ) then N(x) = vr"™/ 2 dn /2 P{x) and hence iV(0 = ir- n / 2 {i £ n )" /2 -P(£)- % 
definition, for all (j) E we have 

<P, 0) = lim <P, X e * 0) = hm (P, (3.1) 

e->0+ e-S-0+ 

Now (all integrals below are absolutely convergent) 

_ /• p -e 2 \x'\ 2 ( \ r poo i- -e 2 \x'\ 2 -ex n . 

(P, = / . dx= / / - G e— " « 0(0 « d*n dx' 

J|x'| 2 <x„ \/Sn — |x'| z Jl"-i J|x'| 2 il" V X n — |x'| z 
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where 



r-oo p -e 2 \x'\ 2 —exn 

Fe(0= I I , ,, 2 e~ lX< dx n dx' 



J 



v {n - 1)/2 ( -in 2 \ 

exp 



(^ e 2 + e + ^ n )n-i ^^4(e 2 + e + ^n)7 y/t+%& 

and the branch of square root used is the one obtained by cutting out the negative reals and whose 
argument lies in (— tt/2, 7r/2). Hence 

= lim / n- n / 2 (iz n r/ 2 F e (om<%- 

Now 7r -n / 2 (i£, n ) n ^ 2 Fe(0 m &y be seen to be bounded with a bound independent of e and its pointwise 
limit is exp(i|£'| 2 /(4$; n )) almost everywhere, hence by the dominated convergence theorem 

{$,</>) = [ exp(i|£'| 2 /(4£„))0(Od£. 

QED 



4 Proof of Proposition [3] 



Because of linearity it is enough to deal only with the two special cases where either / = or g = 0. 
We give the proof only of the case where g = 0, in which case v is given by the formulas (|1.5p . (|1.7p 
for n odd, even respectively. The case where / = is dealt with similarly with the formula for the 
new v being almost the same as (|1.5|) . (|1.7|) except the t term is dropped, the D power is reduced 
by one, and the expressions must be multiplied by 2. 



4.1 Odd n case 



Let / G Co°(IR n ) with / supported in the centered ball of radius R. Since (Mf)(x,t) ^ iff the 
sphere \y — x\ = t intersects the ball \y\ < R, we see that (Mf)(x,t) is supported in the region 
1 1 — \x\ | < R, hence, using (jl.4p . u(x,t) is supported in the region 1 1 — \x\ \ < R. Now suppose 
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t > 0, p a real number and h G Co°(R n ); then 

D(iP(Mh)(x,t)) = ^tP- 2 (Mh)(x,t) + -^— [ 0.{VK){x + te)dO 
2 2\S n l \ J\ e \ =1 

= l t p- 2 (Mh)(x,t) + J^- [ (x + te- x )-(vh)(x + te)de 

2 2 \ s I J\e\=i 

= ci t p ~ 2 {Mh){x, t) + c 2 t p - 2 M(y ■ Vh{y)){x, t) + c 3 ^~ 2 x • M(Vh)(x, t) (4.1) 

= t p - 2 x a M(h a )(x,t) (4.2) 
M<i 

where h a (x) are smooth functions which are just made of derivatives of h, possibly multiplied with 
powers of x. Hence, using (|4.2p (n — l)/2 times in (jl.4j) . we see that 

«(a;,t) = ctD^ n - l ^ 2 t n - 2 {Mf){x,t) = x a (Mh a ){x,t) (4.3) 

M<(n-l)/2 

where h a are made up of derivatives of /, possibly multiplied with powers of x. Now for any h 
supported in |x| < R, the surface area of the intersection of the sphere |y — x\ = t with the ball 
\y\ < R will be bounded above by a constant independent of x or t. Hence for t > 1 

\(Mh)(x,t)\ < I [ \h(y)\dS y 4 + (4.4) 

1 P I J\y-x\=t 

and noting that u is supported in \t — \x\ \ < R, for t > 1 we have from (14. 3D that 

KM)| * (1 ± « (l + |t|)-(^)/ 2 . 



4.2 Even n case 



First note that u(x, t) is supported in \x\ < t + i?. We will show the decay rates, as t — > oo, for the 
two cases \x\ — R < t < \x\ + 3R and |x| + 3R < t separately. 

Suppose \x\ + 3R < t; then for any y in the support of / we have \x — y\ < \x\ + \y\ < \x\+ R < t. 
Hence the region \x — y\ < t includes the support of /. So from (|1 .7j) we have 

M M dy 



i(x,t)\ = ±ctD n / 2 [ fiy) dy = ±ct ! 

J\v\<R \/t z - \x - yr JW 



i\ y \<R V^-k-yl 2 ' J\v\<R (t 2 ~ I* - yl 2 )^ 2 
t t 1 



(t 2 - (\x\ + U)2)(n+l)/2 ( t + | x | + fl)(n+l)/2 (£ _ | x | _ fl)(n+l)/2 ' 

Now for |x| + 3i? < t we have R< (t — \x\)/3 and 

2 

i - |x| - i? > t - |x| - 2i? + i? > t - \x\ - -(t - \x) + R > c(\t - \x\\ + 1), 
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hence 

\u(x,t)\ 4 (l + t)-( n - 1 )/ 2 (l + |t-|x||)- (n+1)/2 . 

This estimate is better than the estimate in the statement of the proposition - the weaker estimate 
in the statement of the proposition is tight when / = and g is non-zero. 

For \x\ — R < t < \x\+ 3R we use a technique similar to one used in the odd case. Since we are 
interested in situation when t is large, we may assume that \x\ > max(l,i?). From (|1.7p we have 

pt n— 1 pi „n— 1 

u{x,t) = ctD n ' 2 (Mf)(x,r)dr = ctD n / 2 t n - 1 / -==(Mf)(x,ta)da. (4.5) 

Jo Vt z — r z Jo vl-« 

Now, for t > 0, for any real number p and any smooth function h{x) we have 

D (t p (Mh)(x, st)) = £ t p ~ 2 (Mh)(x, st) + f n , f sO ■ (Vh)(x + stff) d6 
2 2 |6 n I J\e\=i 

= V - t p ~ 2 (Mh) (x, st) + J** . / (x + - x) ■ (V/i) (s + st0) d# 
2 2|6 n x | 7| | =1 

= ci t p ~ 2 (Mh)(x, st) + c 2 t p - 2 M(y • V%))(x, st) + c 3 t p_2 x • M{Vh)(x, st) 
= t p ~ 2 x a M(h a )(x,st) 

\a\<l 

where h a (x) are smooth functions which are just made of derivatives of h, possibly multiplied with 
powers of x. Hence 

tD n l 2 t n - l (Mf)(x,ts)= ^ x a {Mh a ){x,st); 

\a\<n/2 

using this in (|4.5p . noting that (Mh a )(x, st) is supported in \x\ — R < st < \x\ + R, and using (|4.4|) 
with h replaced by h a , we obtain (note \x\ > R) 

I' 1 s n ~ l \x\ n / 2 f 1 1 

\u(x,t)\ 4 \x\ n ' 2 r \(Mh Q )(x,st)\ds 4 ^rr / n o ds 

JO VI - S 2 t n J(\x\-R)/t VI - s 2 

\ x \n/2 

= -^zr arccos((|x| - R)/t). 

We are considering the region where t — 3R < \x\ < t + R, and there is no loss of generality in 
assuming that t > max(l,4i?); hence 

\u(x,t)\ 4 (1 + t)- (n - 2 ^ 2 arccos((t - AR)/t) 4 (1 + t)" (n " 1)/2 

where the last step follows from lim^oo t 1 / 2 arccos((i — 4R)/t) = c ^ by L'Hopital's rule. So 
Proposition [3] follows from the fact that R < \t — \x\ \ < 3R in the region of interest. QED 
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